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AND OPERATING CHARACTERISTICS OF A BELLOWS -ACTUATED 
SPLIT -FLAP AILERON ON A 60° DELTA WING AT 
MACH NUMBERS BETWEEN 0.8 AND 1.8 
By Eugene D. Schult 


SUMMARY 


A free -flight investigation has been conducted by the Langley 
Pilotless Aircraft Research Division in the Mach number range between 0.8 
and 1.8 to determine the maximum deflection, the zero-lift rolling effec- 
tiveness, and the general operating characteristics of a bellows -actuated 
aileron control system energized by the pitot (impact) pressure. These 
tests of a partial-span split -flap aileron on a 60 ° delta wing indicated 
that this system affords a promising means for obtaining lateral control 
at supersonic speeds in that substantial aileron deflections axe readily 
obtained without the aid of an auxiliary servo power supply. Maximum 
aileron deflections of the order of 20° at a Mach number of 1.8 were 
measured when using a simple one-cell bellows arrangement, and the deflec- 
tion generally increased with increasing supersonic Mach number . Higher 
deflections may be expected with improved bellows design. The maxi mum 
rolling effectiveness of the test configuration remained constant at high 
subsonic speeds; then decreased at supersonic speeds to approximately 
one-third the subsonic level. No aileron flutter was observed. Suitable 
methods for regulating the aileron deflection are described. 


INTRODUCTION 


One of the problems associated with the design of controls for air- 
planes and missiles flying at supersonic speeds is that of providing an 
effective means for overcoming the aerodynamic loads due to control 
deflection. Because of severe space limitations, the use of massive 
control elements or a large power-boost system is not always practical, 
with the result that an increased emphasis has been placed on simplified 
controls which utilize the energy in the airstream to deflect the control 
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surface. By this concept it may be possible to eliminate altogether the 
weight and cost of the servo power supply and storage component. 

A proposed lateral control of this type is the bellows -actuated 
aileron. One possible arrangement (fig. l) consists of a split-flap 
aileron which is deflected by means of a pneumatic bellows, vising air 
obtained from a pitot (impact) pressure source. Early investigations 
with air brakes and landing flaps have shown that large deflections are 
possible with this system (refs. 1 to 4). So far as is known, the first 
mention of applying bellows to lateral or other controls where precise 
control setting is necessary was made in reference which included 
some preliminary data from the present investigation. 

In the present investigation the Langley Pilotless Aircraft Research 
Division conducted a free -flight test of a simple, bellows -actuated ail- 
eron located slightly forward of the trailing edge of a 60° delta wing. 
This aileron arrangement causes little structural interference with 
trailing -edge high-lift devices and provides lateral control at low and 
transonic speeds (refs. 6 and 7); only limited data are available on 
its effectiveness at supersonic speeds, however. Measurements were made 
of the maximum aileron deflection and zero-lift rolling effectiveness at 

Mach numbers between 0.8 and 1.8 and Reynolds numbers between 4 X 10^ 

£ 

and 13 X 10°. Also included are some limited data on aileron hinge 
moments, aileron response time, bellows effectiveness, and the operating 
characteristics of a proposed aileron-deflection control system. 


SYMBOLS 


A 


aspect ratio. 



S 


total wing area, sq ft 


b 


total wing span, ft 


c 


chord, ft 


A 


sweepback angle of wing leading edge, deg 


5 


aileron deflection, deg 


t 


time , sec 


V 


velocity of model, ft/sec 


M 


Mach number 
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R 

P 

2b 

2V 

P 

<1 

P’ 

yf 

m 

“B 


Reynolds number, based on mean aerodynamic chord of wing 
rolling velocity of model, radians/sec 

wing -tip helix angle, radians 

absolute pressure, lb/sq ft 

free-stream dynamic pressure, 0.7P a M^, lb/sq ft 

P - Po 

pressure coefficient, — ^— =• 
aileron span, ft 

moment area of aileron about aileron hinge axis, Cfyf/2, 
cu ft 

moment area of bellows "footprint" on aileron lower surface 
about aileron hinge axis, cu ft 


m B 

m 

H 

c h 

C b6 

N 

D 

l 

d 

r 

w 

F 


moment -area ratio indicating effectiveness or mechanical 
advantage of the bellows 

aileron hinge moment, positive when tending to reduce the 
aileron deflection, ft-lb 

aileron hinge-moment coefficient (per aileron), H/2qm 



d5 


number of bellows cells 

depth of bellows cavity, aileron chords 

length of air duct, aileron chords 

inside diameter of air duct, aileron chords 

radius of control -valve core, aileron chords 

width of control-valve ports, measured along aileron 
hinge axis, aileron chords 

frequency, cycles/sec 
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0 phase -lag angle, deg 

3 aileron amplitude ratio (ratio of aileron deflection at any 

given frequency to aileron deflection at zero frequency) 

k x ,k 2 ,k 5 arbitrary constants 
Subscripts : 

a atmospheric or static (pressure) 

b base (pressure behind a bluff body, referring to a less- 

than-atmospheric pressure sink) 

T pitot (impact pressure) 

B bellows 

U upper surface of aileron 

L lower surface of aileron 

f aileron 


MODEL 


The test vehicle employed in this investigation is illustrated in 
figure 2 and detailed in figure 3. The delta wings had leading edges 
swept back 60° and modified hexagonal sections of constant thickness. 

The ratio of wing thickness to chord varied from 3 percent at the fuse- 
lage juncture to a maximum of 9 percent near the tip. A free-to-roll 
tail assembly stabilized the model longitudinally without introducing 
rolling moments. 

The lateral control system consisted of an essentially free-floating 
split -flap aileron located on the upper surface of each wing and deflected 
by means of single-cell bellows (fig. 3(b)). The bellows were ducted to 
a pitot and base pressure probe through separate control valves (fig. 3(c)). 
A motor-driven cam programed the control valves so that the ailerons were 
deflected independently as functions of time . Control valve (l) was used 
to pulse the pressure in bellows (l) alternately between the inlet and 
base pressures. In addition to deflecting the aileron to obtain roll 
data, this scheme provided an indication of the system time lag over the 
test Mach number range. Valve (2) regulated the deflection of aileron (2) 
in order to provide test data at intermediate deflections . This was 
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accomplished by varying the flow of air through a static-bleed port so 
that the bellows pressure could be adjusted to any desired fraction of 
the inlet pressure. Bench tests of this device indicated that, for a 
given valve setting, the bellows pressure ratio P'b/P't was essen - 
tially a constant within the range of inlet pressures encountered during 
the flight tests . The bellovs pressure ratio did not vary linearly with 
valve setting, however, and in the interest of improving the quality of 
the data it was necessary to modify the cam profile to make the ratio 
vary linearly with cam rotation and time. 

A single probe served as a source of pitot pressure to inflate the 
bellows and as a reservoir of low base pressures to assist in evacuating 
the bellows quickly. 


TEST TECHNIQUE 


The flight test was conducted at the Langley Pilotless Aircraft 
Research Station at Wallops Island, Va. A single-stage booster, con- 
sisting of a 6-inch-diameter ABL Deacon rocket motor, accelerated the 
model at essentially zero lift, yaw, and roll to a maximum Mach num- 
ber of I .85 in 3.0 seconds. During this time the control-valve set- 
tings were fixed in approximately the positions shown in figure 3 (c). 
After the booster separated from the model, the programing motor was 
energized; this caused the ailerons to be pulsed and the model to roll 
alternately to the right and left while decelerating through the test 
Mach number range. 

Measurements were made of the velocity by using continuous -wave 
Doppler radar and of the rolling velocity by using polari zed-wave radio 
equipment. These data, in conjunction with SCR 58 ^- space radar and 
radiosonde measurements, permitted an evaluation of the Mach number M 
and the wing -tip helix angle pb/2V as functions of time. Simultane- 
ous telemeter broadcasts provided time histories of the aileron deflec- 
tions and the bellows pressure, making it possible to determine the 
maximum aileron deflection ^, aY , the rolling effectiveness of the ail- 
eron per degree and the aileron hinge-moment coefficient C^. 

The aileron hinge -moment coefficients, which reflect the hinge moments 
balanced by bellows (2, , were derived by use of the following relation: 


c h 


H _ £j3 ^B 
2qm 2 m 


( 1 ) 
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The moment -area ratio — is essentially the mechanical advantage of 

m 

the bellows and is a constant at a given aileron deflection. Its vari- 
ation with aileron deflection was determined from a bench test of the 
configuration by substituting the experimental quantities of hinge 
moment and bellows pressure into equation (l) and solving for the ratio 
at the observed aileron deflections. Since the bellows was not attached 
to the aileron, it was possible to measure only positive hinge moments 
with this arrangement. 

The test Reynolds number, based on the mean aerodynamic chord of 
the wing, varied between 4 x lO^ and 13 x 10 ^ and is shown plotted 
against Mach number in figure 4 . 

ACCURACY OF MEASURED DATA 


The test results are believed to be accurate within the following 
limits : 


M 

pb/ 2 V, radians 
Ch 

&, deg . . . 


Subsonic 

± 0.01 

±0.003 

±0.03 

± 0.4 


Supersonic 

± 0.01 

± 0.002 

±0.02 

± 0.4 


RESULTS AND DISCUSSION 


The test results of the present investigation are presented in fig- 
ures 5 to 9. In figure 5 flight histories of the measured bellows pres- 
sure, aileron deflections, and model rolling velocity are plotted against 
Mach number. The bellows pressure in cell (2) has been reduced to coef- 
ficient form and presented as a fraction of the calculated pitot-pressure 
coefficient. No pressures were measured in cell (l) during the test. 

The aileron deflection data show that both ailerons floated freely at 
small deflections; this phenomenon, which was especially noticeable at 
high subsonic speeds, is attributed to an overbalanced hinge moment con- 
dition and will be discussed later. No indications of aileron flutter 
were observed over the test Mach number range. 
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Aileron. Rolling Effectiveness 

Figure 6 presents the variations with Mach number of the maximum 
aileron deflection, the rolling effectiveness per degree of aileron 
deflection, and the maximum rolling effectiveness. The results show 
that substantial aileron deflections were obtained with a simple, one- 
cell bellows arrangement and that the maximum deflection generally 
increased with increasing supersonic Mach number. It will be shown 
later that larger deflections may be expected if the number of cells 
is increased. The slightly lower deflection indicated for aileron (T) 
is believed to have been caused by air leakage from valve (l) due to 
construction inaccuracies . 


The rolling -effectiveness parameter 


pb/2V 

6 


was obtained from essen- 


tially steady-state values of rolling velocity measured near the maximum 
deflection of a single aileron. Neglected in this presentation are the 
small losses in rolling effectiveness (0.5 percent) caused by bearing 
friction in the free-to-roll tail. The slight deflections of the oppo- 
site aileron (fig. 5 ) have also been neglected, inasmuch as reference 7 
has indicated that this aileron may be relatively ineffective at small 
angles unless a wing slot and deflector plate arrangement is incorporated. 


At supersonic speeds the rolling effectiveness of the split-flap 
aileron used in this investigation is estimated to be approximately the 
same as that of a conventional sharp -trailing -edge aileron of the same 
span and chord. In this comparison the measured rolling effectiveness 
of a full-span trailing -edge aileron on a 60 ° delta wing (ref. 8) was 
corrected for differences in aileron span and chord by use of refer- 
ence 9* The final plot of figure 6 shows that the maximum rolling effec- 
tiveness of the present-test configuration remained constant at high 
subsonic speeds, primarily because of the variations in the maximum ail- 
eron deflection. At supersonic speeds the maximum rolling effectiveness 
decreased to approximately one -third the subsonic level. 


Aileron Hinge Moments 

Figure 7 presents the aileron hinge -moment coefficients and 

as functions of aileron deflection and Mach number, respectively. These 
coefficients were derived by use of equation (l) from the aerodynamic 
hinge moments neutralized by the bellows. The data for the boosted 
(no-roll) phase of the flight were corrected to essentially steady-state 
load conditions by taking into account the longitudinal accelerations 
acting on the aileron. For the coasting (rolling) phase the total error 
in hinge moment resulting from model deceleration, changes in model 
rolling velocity, and friction in the aileron hinge was estimated to be 
less than 2 percent of the maximum hinge moment at a given Mach number 
and was neglected. A small Mach number correction was made to the data 
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of figure 7(a) to account for the difference between the actual Mach num- 
ber and the average Mach number for the cycle; this correction was gen- 
erally less than 5 percent. The differences between the curves for 
increasing and decreasing aileron deflections are attributed largely to 
the differences in the average rolling velocity for the two cases . 

The test results at high subsonic speeds are in general agreement 
with the data of reference 7 for a similar aileron on a thinner wing 
configuration (t/c = 1.5 to percent). An overbalanced hinge moment 
condition existed at small aileron deflections and caused the ailerons 
to float as shown in figure 5* The condition, which was more pronounced 
at subsonic speeds, coincides with the reversed lift and pitching -moment 
increments observed in other tests of split flaps located forward of the 
wing trailing edge; the phenomenon occurred when the flow reattached to 
the wing surface behind the flap (ref. 10). 

At supersonic speeds the experimental results are in good agreement 
with values calculated from simple theory. The calculated values were 
obtained by use of the following relations: 


c h ~ c hjj + c h L (2) 


c hu - ^ (3) 



The pressure coefficients for the upper surface of the aileron P'y 

were calculated from two-dimensional second-order theory. The base 
pressure coefficient P'^ = -l/M^ corresponds to a limiting value of 

negative pressure behind the aileron (ref. 11) and is assumed to act 
on the lower surface of the aileron not in contact with the bellows. 


System Time Lag 

Since the time required to operate a bellows -actuated aileron may 
be of interest, measurements were made of the time increments required 
to deflect the aileron to approximately 90 percent of maximum steady- 
state aileron deflection and to retract it to approximately zero deflec- 
tion. The time increments for these two cases were essentially the same 
and a curve faired through the test points is presented in figure 8. As 
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a matter of interest, these time intervals are of the same order as those 
generally obtained with conventional control -pulsing mechanisms. 


Estimation of Aileron Deflection 

In appendix A a method is presented for estimating the deflection 
of a bellows -actuated aileron. These estimates are shown in figure 9 
to be in good agreement with experimental values at M = I.76. In fig- 
ure 10 this method is extended to higher Mach numbers in an attempt to 
determine the variation of the maximum aileron deflection with Mach num- 
ber and the additional deflection which could be obtained by adding 
another cell to the bellows. The results indicate that the maximum 
deflection of the aileron used in the present investigation continues 
to increase with increasing Mach number and approaches a limiting value 
of approximately J0° at M = 4.0. Adding another cell increases this 
deflection approximately 30 percent. The base pressure behind the ail- 
eron is estimated to have little effect on the maximum aileron deflec- 
tions at these higher Mach numbers . 


Aileron Deflection Control 

One difficulty encountered with the present method of regulating 
the aileron deflection (by regulating the bellows pressure) was that 
of maintaining precise control over the deflection. As a result, a 
system incorporating feedback was devised and preliminary bench tests 
were conducted (appendix B) . The results of these tests were encouraging 
and the system will be used in later research with a pitch control 
employing a bellows -actuated flap. 

CONCLUSIONS 


A free -flight investigation was conducted by the Pilotless Aircraft 
Research Division in the Mach number range between 0.8 and 1.8 to deter- 
mine the maximum aileron deflection, the zero-lift rolling effectiveness, 
and the general control characteristics of a bellows -actuated aileron 
control system energized by the free-stream pitot pressure. The fol- 
lowing conclusions are based on tests of an inboard, partial-span, split- 
flap aileron located slightly forward of the trailing edge of a 60° delta 
wing: 

1. Maximum aileron deflections of the order of 20° at M = 1.8 were 
measured with a simple one -cell bellows, and the deflection generally 
increased with increasing supersonic Mach number. Calculations show that 
the deflection may be increased further by slight improvements in the 
bellows design. 
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2. The maximum rolling effectiveness remained constant at high sub- 
sonic speeds and then decreased at supersonic speeds to approximately 
one-third the subsonic level. At supersonic speeds the effectiveness of 
this aileron was estimated to be the same as that of a conventional 
sharp-trailing-edge aileron of the same span and chord. 

3 . Slightly overbalanced hinge moments were observed at low deflec- 
tions for this aileron configuration. The phenomenon was more pronounced 
at subsonic than at supersonic speeds . No indications of aileron flutter 
were observed over the test Mach number range . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., July 29 , 195^- • 
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APPENDIX A 

A METHOD FOR ESTIMATING AILERON DEFLECTION 
AT SUPERSONIC SPEEDS 


A steady-state deflection of the bellows -actuated aileron occurs 
when the hinge moment available for balancing (due to the bellows) is 
equal to the hinge moment required. The hinge moment required is, under 
static load conditions, essentially equal to the sum of the hinge moments 
contributed by the aerodynamic pressures acting on both the upper surface 
of the aileron and the lower surface not in contact with the bellows; 
that is, 


C h B " % + 


c h T 


(Al) 


where 



p, B 
2 m 


(A2) 


c hy - 



(A3) 


C *L 




(A4) 


m-Q 

The moment-area ratio or bellows effectiveness parameter — is a func- 
tion of the aileron deflection and can be either estimated or determined 
experimentally from bench tests of a given bellows configuration. An 
estimation of the ratio for cell -type bellows acting on ailerons with 
rectangular plan forms is given by the equation 


m B 

m 


L-j 


(A5) 


where 


k l " 




(A 6) 
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The term kq is, for N cells, the calculated moment-area ratio for an 
infinitely thin, nonstretching, flexible bellows maintained by internal 
pressure in contact with an infinitely long aileron. The corrections 
to account for the losses in bellows effectiveness associated with 
finite bellows thicknesses (cavity depth D) and length-to-chord 


ratios Af are ] 
corrections are : 


k 2 = 


terms 

k 2 

and ] 

_ 1 - JB- 

A + 

n5 V 

2N 


i 8 onJ 


respectively. These 


n 


(A7) 


kx = 1 — 


4 

3A f 


(l 


k l k 2 


(Af ^ 2) (A8) 


The following illustration shows the assumed bellows "footprint" used in 
deriving equations (A6) to (A8). 




mt> 

A comparison of the calculated ratio with values experimentally 

derived for the present-test configuration is shown in figure ll(a). 

The experimental values were slightly lower at the higher aileron 
deflections. Figure 11(b) illustrates the fact that the effectiveness 
of the bellows may be improved significantly by increasing the number 
of cells. 

Using the calculated moment-area ratio (eq. (A5)), the aileron 
deflections of the test configuration were estimated for M = 1.76 by 
means of equations (Al) to (A4). The pressure coefficient P'y was 

calculated from two-dimensional second-order theory, and P'q, the lim- 
iting negative pressure coefficient behind the aileron (p'^ = -1 /m^) , 

was obtained from reference 11. The results are presented in figure 9 
and show good agreement with experimental values. 
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APPENDIX B 

AILERON DEFLECTION CONTROL 

In the test configuration the aileron deflection was regulated by 
proportioning the pressure in the bellows between the static and the 
pitot pressures (fig. 5)» This simple method is generally not suitable 
for maintaining precise control over the aileron, especially when the 
maximum deflection varies with Mach number and when nonlinear hinge 
moments are likely to occur. One method of obtaining precise aileron 
deflection control is to provide the system with feedback (fig. 12), 
so that any variations from the desired deflection, which may result 
from an increase or decrease in hinge moments, will be compensated for 
by the maximum positive or negative pressures available. 

A mock-up of a bellows -actuated aileron based on the system shown in 
figure 12(a) is illustrated in figure 13, and results of static response 
tests are presented in figure lk. The aileron hinge-moment loading was 
simulated for M = 1.9 and the inlet pressure used was the pitot pres- 
sure at this Mach number. The results show that, with full pitot pres- 
sure available ^P'/P'fj = 1.6) , the output is linear with input until the 
maximum aileron deflection at this Mach number is approached (21°). The 
difference between the experimental and perfect response curves was caused 
largely by air leakage resulting from valve construction inaccuracies. A 
good indication of the control "stiffness" may be gained by comparing the 
curves for P'/P'rp = 1.0 and P'/P't = 0*7> here a 30-percent reduction 
in the available pressure for the same aileron loading is analogous to 
approximately a kO-percent increase in aileron hinge moment when 
P'/P't = 1*0. Similarly, a 70 - P erc ent reduction in the available pressure 
for the same loading is analogous to a 300-percent increase in hinge 
moment . 

Some limited data on the frequency-response characteristics of this 
configuration are presented in figure 15. These tests were of an explora- 
tory nature, not intended to determine an optimum configuration, but to 
illustrate the effects of variations in several design parameters on the 
control response . Due to the method of applying the load, some friction 
damping was introduced which may have altered the general shape of the 
curves somewhat; this should not, however, impair the usefulness of the 
data. In the presentation the amplitude ratio is defined as the ratio 
of the aileron deflection at any given frequency to the aileron deflec- 
tion at zero frequency. 

The results show that an increase in input amplitude from 6° to l6° 
(reflecting a 60-percent increase in bellows volume) or changes in the 
available pressure or duct dimensions produced only slight changes in 
the response of the configuration. The greatest improvement in the fre- 
quency response was obtained by increasing the valve port width w. 
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Figure 1.- Schematic sketch of a possible wing -aileron arrangement in 
which bellows are employed to actuate the ailerons . 
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(a) Model configuration. 
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Figure 2.- Photographs of configuration used in the present investigation. 
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(b) Wing configuration disassembled to show control mechanism. 
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Figure 2.- Continued. 
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(c) Model and booster on launching stand. 


Figure 2.- Concluded. 
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Free-to-roll tail 
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Figure 3 .- Geometric details of configuration used in the present inves- 
tigation. All dimensions are in inches. 


NACA RM L54H17 CONFIDENTIAL 


CONFIDENTIAL 



(b) Wing and aileron arrangement . 
Figure 3*~ Continued. 
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Boosted phase 



M 

Figure 4.- Variation of test Reynolds number with Mach number; Reynolds 
number is based on wing mean aerodynamic chord ( 1.158 ft). 
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Boosted phase (no-roll) 

Coasting phase 
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Figure 5*- A flight history, plotted as a function of Mach number, of 
measured values of bellows pressure in cell (2), aileron deflec- 
tions, and resulting rolling velocities. 
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Figure 6 . - Variations with Mach number of measured maximum aileron 
deflection and rolling effectiveness obtained with the bellows - 
actuated aileron control configuration of the present test, and 
comparison with other data. 
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(a) Typical variations with aileron deflection of the hinge-moment coef- 
ficient neutralized by the bellows, and comparisons with experimental 
data from reference 7 and two-dimensional second-order theory. 



M 
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(b) Variation with Mach number of average slopes of hinge -moment- 
coefficient curves for the present-test aileron. 

Figure 7*- Measured hinge-moment characteristics of the test aileron. 
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Figure 8.- Variation with Mach number of maximum measured time to deflect 
or retract aileron and of chord lengths traveled during this interval. 
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(a) Calculation of the maximum aileron deflection obtained with pitot 

pressure . 



(b) Calculated aileron deflections for bellows pressures less than pitot 

pressure. 

Figure 9.- Comparison of calculated and experimental aileron deflections 
obtained with bellows pressures up to pitot pressure M = I.76; 

Af = 4; D = 0.1. 
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Figure 10.- Variation with Mach number of calculated maximum aileron 
deflection for a single -cell and a double -cell bellows inflated 
with pitot pressure Af = 4.0j D = 0.1. 
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(a) Variation with aileron deflection of calculated and experimentally 
derived values of bellows effectiveness parameter N = 1. 



6, deg 


(b) Effect of number of bellows cells on calculated bellows effective- 
ness parameter A^ = 4; D = 0.1. 

Figure 11.- Variations with aileron deflection of the bellows moment-area 

ratio or effectiveness parameter . 
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(a) Control employing rotary-type valve. 


Input 



(b ) Control employing lineal- type valve. 

Figure 12.- Two variations of an aileron-deflection control employing 

feedback. 
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Figure 13*- Mock-up of bellows -actuated aileron with deflection control L- 8 II 56 
employing mechanical feedback. Aileron plan form and bellows are 
identical to those in flight-test model. 
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Figure 14.- Static response characteristics of an aileron-deflection 
control with feedback. Aileron loading H/6 and pitot pressure Prp 

are simulated for M = 1.9 at sea level, r = 0.1; w = 0.15; Ap = 4; 
D = 0.1. 
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(a) Same effects of variations in input amplitude & input and inlet 
pressure ratio P'/P'ij on the phase-lag angle and output ampli- 
tude ratio- d = 0.1; 1=4; w = 0 . 15 . 



(b) Some effects of changes in duct and valve-port geometry on the 
phase -lag angle and output amplitude ratio 6 input = 10 °j 
P'/P' T = 1.0. 

Figure 15*- Frequency-response characteristics of a bellows -actuated 
aileron-deflection control with feedback. Aileron loading and 
pitot pressure simulated for M = 1.9 at sea level, r = 0.1; 

Af = 4; D = 0.1. 
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